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Recent theories of Sr 2 Ru 04 based on the interplay of strong interactions, spin-orbit coupling and multi-band 
anisotropy predict chiral or helical ground states with strong anisotropy of the pairing states, with deep minima 
in the excitation gap, as well as strong phase anisotropy for the chiral ground state. We develop time-dependent 
mean field theory to calculate the Bosonic spectrum for the class of 2D chiral superconductors spanning 3 He-A 
to chiral superconductors with strong anisotropy. Chiral superconductors support a pair of massive Bosonic ex¬ 
citations of the time-reversed pairs labeled by their parity under charge conjugation. These modes are degenerate 
for 2D 3 He-A. Crystal field anisotropy lifts the degeneracy. Strong anisotropy also leads to low-lying Fermions, 
and thus to channels for the decay of the Bosonic modes. Selection rules and phase space considerations lead 
to large asymmetries in the lifetimes and hybridization of the Bosonic modes with the continuum of un-bound 
Fermion pairs. We also highlight results for the excitation of the Bosonic modes by microwave radiation that 
provide clear signatures of the Bosonic modes of an anisotropic chiral ground state. 


Introduction 

Superfluid 3 He and unconventional superconductors share a common and fundamental property that the ground state breaks 
one or more symmetries of the normal Fermionic vacuum in conjunction with the usual U(l) gauge symmetry associated with 
BCS condensation. In the case of Sr 2 Ru 04 the connection with superfluid 3 He may be stronger. The theoretical proposal for 
a spin-triplet, p-wave ground state in Sr^RuC^ was motivated by similarities between Sr 2 Ru 04 and liquid 3 He, particularly the 
existence of exchange enhanced paramagnetism in a strongly correlated Fermi liquid fT]. In liquid 3 He long-lived ferromagnetic 
spin fluctuations (“paramagnons”) are believed to be the mechanism responsible for the BCS pairing instability to a spin-triplet, 
p-wave ground state (2] [3). The Balian-Werthamer (BW) state (4), identified as the B-phase, with total angular momentum 
/ = 0 fully gaps the 3D Fermi surface, and as a result minimizes the free energy in the weak-coupling limit. The BW state 
is time-reversal invariant, but spontaneously breaks relative spin-orbit rotation symmetry. However, spin-triplet correlations, 
which differentiate between the Anderson-Morel (AM) state and BW states, feedback to modify the spin-fluctuation exchange 
interaction, leading to stabilization of the AM state at high pressures mm- This is the chiral A-phase which breaks time-reversal 
(T) symmetry and reflection symmetry in any plane containing the chiral axis (P 2 ), but preserves T x P 2 (chiral symmetry). 

Similarly, Rice and Sigrist argued that for S^RuCH, which is a layered perovskite belonging to the D 4 h point group, the likely 
pairing state, in analogy with 3 He, would be the planar state, the 2D analog of the BW state. In weak-coupling theory, and ne¬ 
glecting spin-orbit coupling, the planar state, which belongs to the one-dimensional Aj u representation of D 4 h, is degenerate with 
the 2D chiral AM state, which belongs to the 2D E| Ll representation. Thus, if spin-fluctuation exchange is also the mechanism 
for pairing in S^RuCU, we expect the strong-coupling feedback effect will stabilize the chiral AM state. However, in addition to 
the near two-dimensionality of the Fermi surface of SnRufTt, incommensurate spin-density-wave (SDW) fluctuations connected 
with the quasi-one-dimensional a and /3 bands co-exist with ferromagnetic fluctuations at lower temperatures 0- Spin-orbit 
coupling, and the possibility of pairing on multiple Fermi surface sheets likely play important roles in determining the pairing 
symmetry class, ground state order parameter ifST- TiD as well as the Bosonic excitation spectrum in S^RuCV 

In what follows we develop a field theory for the Bosonic modes based on the Ginzburg-Landau theory (TDGL) for spin- 
triplet, odd-parity pairing that allows one to compare predictions for 3 He and Sr 2 Ru 04 , and to examine the roles of crystalline 
anisotropy and strong correlation effects on the Bosonic modes. The Bosonic modes are excitations of a condensate of Cooper 
pairs for which the parent state is the Fermionic vacuum. Thus, although the TDGL field theory provides insight into the Bosonic 
excitation spectrum, it misses key features of a more complete theory of the low-lying excitations of the BCS pair condensate. 
Notably, (i) polarization effects of the underlying Fermionic vacuum on the excitation energies (masses and dispersion) of the 
Bosonic modes, (ii) finite lifetimes of the Bosonic excitations due to coupling to the Fermionic continuum and (iii) selection rules 
and matrix elements for the coupling of the Bosonic modes to the electromagnetic (EM) field. Key predictions of a microscopic 
theory of interacting Fermionic and Bosonics modes are summarized, including microwave spectroscopic signatures of the 
Bosonic excitation spectrum. 
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Order Parameter 

The mean-field order parameter for spin-triplet Cooper pairs, A a p (p) = d u (p)(io^O y ) a p, where oc.fi label the projections of 
fermion spins of the Cooper pair, {/ay, O y \ jl = x' ,/, z 7 } is the spin-triplet basis of 2 x 2 matrices, is expressed as 

d^i( p)= Y A^ipi), /I £{x',y',z 1 }, ( 1 ) 

i=x,y,z 

where t/ u (p) is a vector under rotations in spin space, and is a function of the vector basis of p-wave orbital basis functions, 
(Pc,p-y, p.), for bulk 3 He. Thus, the amplitudes 4^, provide a bi-vector representation of S0(3) s x SQ(3) L . For quasi-two- 
dimensional Sr 2 Ru0 4 , the orbital basis, (K v (p). K v (p)), provides a 2D vector representation in which Y XJ { p) transforms as p Yy 
under the point group D^. These basis functions reflect anisotropic pairing on the Fermi surface. 

The chiral AM state is the stable equilibrium phase in a narrow temperature window 7 ab < T <T C of bulk ’He at pressures 
P>P C = 21 bar. However, if we confine ’He as a thin film or within a thin cavity of thickness I) < D, ~ 1 Li m the quasi-2D 
A-phase with clip) = -4= % (p A ± ip y ) is the ground state for all pressures and temperatures [; 13]]. We refer to this order parameter 
as the ‘Anderson-Morel (AM) state’, the ‘chiral state’ or the ‘A-phase order parameter’. 

However, in the weak-coupling limit the chiral AM phase is degenerate with the 2D planar phase, d{ p) = (x p t + y p v ). 

In the context of Sr 2 Ru0 4 the 2D planar state is referred to as the “helical state” DU. The degeneracy between the chiral 
and helical ground states is lifted by strong-coupling corrections to weak-coupling theory, or spin-orbit coupling. For 3 He 
strong-coupling effects dominate the nuclear spin-orbit coupling, with the spin-triplet feedback effect on the ferromagnetic spin- 
fluctuation exchange interaction favoring the chiral AM state over the helical state. Recent NMR experiments on ’He confined 
in thin slabs imply that the chiral AM state is stable relative to the helical state down to P ~ Obar m. If superconductivity 
in Sr 2 Ru0 4 is also driven by ferromagnetic spin-fluctuations then we expect the chiral phase to be favored. However, recent 
theoretical calculations based on RG analysis that includes spin-orbit coupling within a multi-band pairing model also find both 
helical and chiral phases, depending upon the interaction parameters defining the lattice pairing Hamiltonian DU- For the chiral 
state the pairing gap, |A(p)| = c/(p)|, is found to be strongly anisotropic on all three bands with deep gap minima over a wide 
range of material parameters DU- Thus, the main observations are: (i) there are competing low energy scales that determine 
the ground state for Sr 2 Ru0 4 , (ii) it is not currently settled whether or not Sr 2 Ru0 4 is a chiral superconductor fl5l . or even 
if the order parameter for Sr 2 Ru0 4 belongs to a two-dimensional representation liT6l . and (iii) whatever the ground state - e.g. 
chiral or helical - low-energy Bosonic excitations may provide unique signatures of the ground state based on their symmetry 
and selection rules. 


Ginzburg-Landau Theory 

Consider the Ginzburg-Landau theory for the class of equal-spin pairing (ESP) states of confined superfluid 3 He in 2D. A 
simple generalization also describes spin-triplet superconductivity on a 2D cylindrical Fermi surface, i.e. pairing on the y band 
in Sr 2 Ru0 4 within the Ej u representation of D 4 | v The general form of the order parameter is then given by 

d (p) = d (A. v Y x (p) +A y Y y (p)) , (2) 

where { Y/(p)\i =x,y} are basis functions for the 2D irreducible representation, E] u , of D 4 | 1 , and A = A v x-t-A y y is a complex 
vector describing pairing in this generalized “p-wave” orbital basis. In what follows we consider the class of ESP states in which 
the direction d along which the Cooper pairs have zero spin projection is fixed, either as a spontaneously broken symmetry 
direction, or by spin-orbit coupling. The general form of the GL free energy functional is constructed from invariants of the 
maximal symmetry group from products of A, and A* through 4th order 03, 

^[A] = Jy dV {a(r) |A | 2 + fa |A | 4 +fi 2 |A • A| 2 + ft [|A . V | 4 + |A_ V | 4 ] } , (3) 

where tt(T) determines the onset of pairing in the Ei u representation, i.e. a ( T,) = 0, and thus, a(T ) « a'(T — T,). The 
fourth order coefficients, j 8 i j 2 j 3 , determine the magnitude of the condensation energy density and the relative stability of phases 
belonging to the Ei u symmetry class. The first three terms in Eq. |3j are invariant under the larger group, i.e. U(l)x S0(2) L; x 
Z“' blt x P x T, while the last term is an additional invariant that is allowed in the case of Sr 2 Ru0 4 by the lower symmetry, 
U(l)x D 4h . 

It is convenient to express the bulk order parameter in terms of an amplitude and a normalized complex vector, A = A a, 
with |a | 2 = a ■ a* = 1. The normalized order parameter is then parametrized by two angles, a = (cos Px + e' <f> sin Av) / y/2 The 
resulting GL functional is then 

^■[A,a] = dV {a(T)A 2 + j3[a] A 4 } , with j3[a] = J3j + /3 2 |a• a | 2 + /3 3 [\a x \ A + |a v | 4 ]. 


(4) 




FIG. 1: Left: Ginzburg-Landau phase diagram for Ei u pairing. The weak-coupling BCS theory predicts the jS parameters lie on the dotted 
red line depending on the degree of anisotropy of the Ej u basis functions, implying the chiral AM state is stable. Right: Anisotropy of the GL 
P parameters in the weak-coupling limit based on the anisotropic pairing model defined in Eq. [l4]as a function of the anisotropy parameter e 
calculated from Eqs. eh m 


For T < T c , minimization with respect to A leads to A^ in = —a(T)/ 2/3 [a], and a condensation energy given by [a] = 
—a(T) 2 / 4/3 [a], with j3[a] >0 for global stability. The ground state is then determined by the normalized order parameter 
that minimizes /3 [a] = j3i + /3i(l — sin 2 (p sin 2 (2$)) -(-/ 33 (1 — ^ sin 2 (2$)). In the case of 2D 3 He we have /I 3 = 0, in which case 
there are two possible ESP ground states; for /3i > 0 and ^2 > 0 the chiral AM state which breaks time-reversal and 2D parity 
is preferred. There are two degenerate chiral ground states which are time-reversed partners of one another, a± = (x ± iy) / y/2. 
However, for — /3i < /T < 0 the in-plane polar state with a = cos fix + sin fly,for. i) G [0, 2k] minimizes the GL functional. 
This phase has a continuous degeneracy with respect to orientation of polar axis in the x — y plane. Tetragonal anisotropy lifts 
the continuous degeneracy of the in-plane polar state. For /T < 0 and > 0 the polar state aligns along a [110] direction, e.g. 
a = (x + y)/\/2. However, for —1(1 + P 2 /P 1 ) < P 3 /P 1 < 0 the polar state aligns along a [100] direction, e.g. a = x. The phase 
diagram is shown in Fig. [T| Note that the weak-coupling prediction for the /j parameters favor a chiral ground state independent 
of the measure of anisotropy, i.e. fo. Furthermore, substantial corrections to weak-coupling theory are required to stabilize 
an in-plane polar state. In fact the helical state, which belongs to a different irreducible representation (A] u ), is a more likely 
competitor to the chiral state since these two states are degenerate in the weak-coupling limit without spin-orbit coupling. In 
addition, AFM spin-fluctuations can lead to attractive, sub-dominant, pairing interactions in even parity, “d-wave”, channels. 
Here we consider fluctuations within the Ej u representation, and neglect possible low-lying fluctuations of the “helical” phase 
(Aj u ) or even-parity, Bj g or B 2 g , “d-wave” states. 


Time Dependent GL Theory - Fluctuations 

Consider the Bosonic excitations of the chiral AM ground state with a+ = (x + iy)/y/2. Space-time fluctuations of the Ei u 
order parameter, = A — Aa , are represented by two complex amplitudes, 

sf{r,t) =D(r,t)a + +E(r,t)a_, (5) 

There are two classes of excitations - modes with chirality L z = +1 represented by the field D(r,t) and modes with the time- 
reversed chirality, L z = — 1, represented by the field E(rj) - and altogether four orbital collective modes within the Ei u represen¬ 
tation. We construct an effective Lagrangian by expanding the GL free energy functional about the ground state. Time-dependent 
fluctuations introduce an additional invariant, GXf = f v dV \ll v/ r G/*, where 11 is the effective inertia of the Cooper pair fluctua¬ 
tions and .(-/ = r) r (/. The effective potential is obtained by expanding the GL functional to 2 nd order in the fluctuations .e/ (r.t): 
'i/ [■<-/] = ■¥\lx \ — .ST[Aa ._]. The Lagrangian, ££ — JC — T/, takes a simplified form when expressed in amplitudes for the normal 
modes: D ± =D±D* and E ± = E ±E*, 


L“ v {g‘ 


( D +) 2 + ( D -) 2 + (£+) 2 + (£-) 2 


— 4A 2 


& (D + y+P2 ((E + y +(£-) 2 )+-MD+r+(£ + n 




( 6 ) 
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Mode 

Symmetry 

Mass 

Name 

EM active 

D+ 

L, = +lC = +l 

2A 

Amplitude 

no 

D 

L, = +lC=—1 

0 

Phase Mode 

no 

E+ 

L- = —1C = +1 

V2A 

E + Mode 

yes 

E~ 

L, = -lC=-l 

y/2A 

E~ Mode 

yes 


TABLE I: Bosonic Mode Spectrum for the 2D Chiral Ground State a+. The masses of the modes are those for an isotropic 2D chiral 
condensate in the weak-coupling limit. Also indicated is whether or not the mode can be excited by absorption of microwave photons. 


The Euler-Lagrange equations reduce to four uncoupled mode equations, 

D c +M; C D C = 0 and E C +M 2 C E C = 0, 


(7) 


where Ml, c is the excitation gap (“mass”) of the Bosonic mode with quantum numbers L z . C, and C is the parity under charge 
conjugation (particle-hole) symmetry mm. The amplitude D , which is a cyclic coordinate in the Lagrangian, is a Goldstone 
mode associated with the broken U( 1) gauge symmetry. Indeed for small amplitude and phase fluctuations of the chiral ground 
state, A(r,f) = A(r,f)e ,e ( r, 0a+ ~ Aa+(1 + (5(r,f) + «0(r,f)), we identify D~ = iA6(r,t) with the phase fluctuation and D + = 
A5(r,l) with the amplitude fluctuation. The amplitude mode of the ground-state order parameter is the Higgs mode for the 
chiral ground state. In particular, D + corresponds to an excitation with the same quantum numbers (L z = + 1 and C = +1) as 
those of the ground state. In the weak-coupling limit the mass of the amplitude mode is determined the pair-breaking energy 
for dissociation into two Fermions, i.e. M++ = 2nip where mp = A is the gap (mass) in the Fermionic spectrum of the broken 
symmetry phase mm. Furthermore, since the amplitude mode has the same quantum numbers as the condensate, the mass of 
the amplitude mode is unshifted, relative to that of two dissociated Fermions, by polarization effects of the underlying Fermionic 
vacuum. We use this to fix the inertia term in the effective Lagrangian as jx = (2/3i 4 - fo), and thus the mass scale for all the 
Bosonic modes of the effective Lagrangian. 

Collective modes of the Cooper pair condensate with quantum numbers distinct from those of the ground state were discussed 
soon after the formulation of BCS theory by Anderson 1221 . Bogoliubov, Shirkov and Tolmachev fl23H . Tsuneto l24l . Vaks, 
Galitskii and Larkin ||25l , Bardasis and Schrieffer [26], and Vdovin [27l . Generalizations of the amplitude mode were discovered 
theoretically in the context of superfluid He by Maki l28l . Wolfle [29], Sauls and Serene [30l . The observation of these 
Bosonic modes using acoustic spectroscopy played an important role in identifying the symmetries of the superfluid phases of 
3 He ED [32]. In this context the collective modes for E| u pairing symmetry with the time-reversed chirality (L z = — 1), both of 
which are massive, correspond to the “clapping modes” of superfluid ’Hc-A in the weak-coupling limit of an isotropic 2D chiral 
AM state (See Table III in Ref. [291 ). In particular, the masses of the E ± modes in the effective Lagrangian are given by 


M_ 


= 2A 


/ + ft 

\2j3i + pi 


1 

2 


M. _ = 2A 


2fe 

2/3 1 +p3 


i 

z 


h=0 

Shlfh=\ 

_fe=o_ 

fh./Pi=k 


> V 2 A, 

(E + mode), 

( 8 ) 

■>V 2 A, 

(E~ mode). 

(9) 


Note that in the weak-coupling limit (P 2 /P 1 = 5 ) for an isotropic (p:\ = 0) 2 D chiral ground state the time-reversed modes are 
degenerate, and lie well below the Fermionic continuum edge at 2A. However, the degeneracy of the E ± modes is lifted by 
tetragonal anisotropy of the Fermi surface and pairing basis functions, which leads to fh / 0. The crystal field splitting of the 
E ± modes is shown in Fig. [ 2 ] In the left panel the splitting of the modes in the weak-coupling limit (P 2 /P 1 = \) is plotted as a 
function of Pi/Pi- The soft E 1 mode, i.e. M. + — > 0 for p\ —> —2pi, is the dynamical signature of the boundary of the unstable 
region of the GL phase diagram shown in Fig. |T| 


Weak-coupling GL Theory for Anisotropic E| u pairing 

For quantitative predictions of the effects of anisotropy of the Fermi surface and pairing interaction on the collective mode 
spectrum we require the angle-resolved density of states on the Fermi surface, n( p), and the anisotropy of the pairing basis 
functions, {7 x (p),7 v (p)}, that define the Ei u orbital order parameter in momentum space, A(p) = A X Y X { p) + A } 7 v (p). For a 
single-band Fermi surface the GL functional for ESP states in the weak-coupling limit can be expressed in terms of Fermi- 
surface averages of the mean-field order parameter defined on the Fermi surface I33ll34l . 

^wc = / v£ /v{a(7')(|A(p)| 2 )+/3o(|A(p)| 4 )} ) 


(10) 
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FIG. 2: Left: Masses of the E± modes versus the GL anisotropy parameter, in the weak-coupling limit, P 2 /P 1 = 5 

based on the weak-coupling j3 parameters defined in Eqs. 1 1| 13 for the anisotropic Ej u basis functions defined in Eq. |14| 
minimum of the pair-breaking edge, 2|A(p)|, bound the gray shaded region. 


. Right: E± masses 
The maximum and 


where (x(T) = In (T/T c ), A> = 7£(3)Nf / (4nT c ) 2 , Nf is the single-spin density of states at the Fermi energy, and (...) = 
f d 2 pn{ p)(...) is the angle-average over the Fermi surface. The pairing basis functions are normalized with respect to the 
anisotropic Fermi surface, (Yf (p)(pj) = 8 , j for i. j £ {x. v}. Evaluating the angular averages in Eq. 
results for the GL material parameters in the weak-coupling limit. 
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gives the following 


pr = 2A)( 


pr = 
pr = 


2A)(|Fv(p)| 2 |f v (p)| 2 
Po{\Y x (p)\ 2 \Y y (p)\ 2 ), 
/3 0 (|F ;c (p)| 4 -3|^(p)| 2 |F y (p)| 2 


( 11 ) 

( 12 ) 

(13) 


These results are obtained by using the transformation properties of the basis functions under the group operations of D^: 

specifically, Y x —> Y y , Y y --T — Y x , Y x —> —Y x , Y y ——> Y y . A key result is that the ratio of GL parameters, pf c /pf c = ], is 
independent of anisotropy 03- Thus, the chiral AM ground state with broken time-reversal symmetry is favored independent 
of the anisotropy of the Ei u basis functions and Fermi surface anisotropy in the weak-coupling limit. 


Cylindrical Symmetry 

For cylindrical symmetry, e.g. a Fermi disk for 2D 3 He-A, we have n(p) = and d 2 p = pj d (!). The normalized p- 

wave basis functions are then Y x = \/2p x = y/2cos0 and Y y = \/2p v = v^sin @, where (j) is the azimuthal angle defining 
the unit vector p normal to the edge of the Fermi disk. Thus, for the chiral ground state the excitation gap (Fermion mass), 
|A(p)| = f |T x (p) —|— iY y (p)| = A, is also isotropic. These basis functions also lead to the following results for the Ei u GL j3 
parameters, A = /3o, P 2 = ' Ai, and ft] = 0 , the latter as expected for an isotropic Fermi surface with pure p-wave basis functions. 
These values give the weak-coupling, isotropic results for the E* modes reported in Tab I ell] i.e. the degenerate “clapping” modes 
of 2D 3 He-A with AT + = M. r = V2A. Note that the Nambu-Sum rule, c = (2A) , is obeyed for both classes (L z = ±1) 

of Bosonic modes for 2D 3 He-A in the weak-coupling limit f35l . 
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Phase Anisotropy 

Recent calculations of pairing instabilities for odd-parity pairing in SriRuCL starting from lattice models based d-band elec¬ 
trons and holes with on-site and near neighbor Hubbard interactions and spin-orbit coupling lead to a helical or chiral ground 
state, but with strong anisotropy. For the chiral state obtained from pairing on hybridized, quasi-ID a and /3 bands (9J, the phase 
of A(p) changes rapidly upon crossing the [110] planes, leading to a low-energy collective mode H36| . To illustrate the effects of 
strong phase anisotropy on the collective modes consider an extreme limit in which the Ei u basis functions are constant within 
each quadrant of the Fermi surface, but change sign abruptly upon crossing any of the [110] planes, i.e. Y x = sgn(p v — p v ) and 
Y y = sgn(p v + p v ). These basis functions have constant amplitude, but are discontinuous in phase (±x) across the [110] planes. 
They give the following results for the Ei u GL /3 parameters, f5\ = 2/3o, /T = /3o, and = — 2/3o. For the E ± modes this leads 
M. . — > 2A and M. + — > 0. The soft E mode is a dynamical signature of the approach to the unstable region of the GL phase 
diagram for J 83 —> — 2/32 as shown in Fig. [I] 


Amplitude Anisotropy 

Multi-band models leading to a chiral superconducting state also exhibit strong anisotropy of the amplitude of the order 
parameter, in particular, an excitation gap |A(p)| with nodal, or near nodal, directions on the Fermi surface lITOl fl~H. Strong 
anisotropy of the Cooper pair amplitudes in momentum space leads to physical properties that are quite distinct from the pre¬ 
dictions based on cylindrical symmetry of the y hand of SroRuO^ including the splitting of the £ ,± modes, and in some cases a 
low-energy collective mode. To illustrate the effects of anisotropy of the pairing interaction on the collective mode spectrum, as 
well as signatures of a chiral ground state with strong anisotropy, we consider the following model for the anisotropic Ej u basis 
functions defined on the y band, or hybridized a and /3 bands, 

Y XJ { p) = V 2 p x , y I(p ), with 7(p) = (1 + e(| 2 p. v p Y | - 1 )) / (1 + 4e(l - e)/n - 2 e(l - 3e/4 ))2 , 0 < e < 1 . (14) 


Tetragonal anisotropy is parametrized by the function /(p), which is invariant under D,^. The Ei u basis functions are normalized, 
/ ^ J / i *(p)T / (p) = Sjj, and in the limit e —>■ 1 exhibit strong anisotropy with a minimum excitation gap, A mm °c (1 — £) —> 0 as 
e —> 1 in the [100] directions. The anisotropy of the excitation gap for the chiral ground state is shown in the left panel of Fig. 
[3] while the corresponding /3 parameters calculated from Eqs. 11 - 13 are shown in the right panel of Fig. [T] The amplitude 
anisotropy has no effect on the ratio, that determines the stability of the chiral state, but has a strong effect on the ratio, 

P 3 /P 1 , that determines the effective potential for the E + modes, and thus the splitting of the these modes by anisotropy, as 
shown in the right panel of Fig. [ 2 ] Note that in addition to the splitting of the E + modes the masses of the E + modes cross the 
continuum edge (2A m j n ) of unbound Fermion pairs. Thus, we expect the E ± modes to become resonances with finite lifetimes 
for sufficiently strong gap anisotropy. However, the theory for the lifetimes of the £ ,± modes is outside the TDGL Lagrangian 
for the Bosonic spectrum, and requires a microscopic theory of the correlated Fermionic vacuum, including the mechanism 
and effects of spontaneous symmetry breaking, and most importantly the interaction and coupling between the Fermionic and 
Bosonic excitations of the chiral superconducting phase. 


Beyond TDGL 

A microscopic formulation of the electrodynamics of the excitations of p-wave superconductors, including the coupling of 
Bosonic modes to a transverse (EM) electromagnetic field, is developed in Refs. f37l[38l for 3D isotropic Fermi systems with 
p-wave, spin-triplet pairing. For a 3D chiral p-wave superconductor in the clean limit the coupled set of linearized dynamical 
equations for the Bosonic mode spectra, including the reaction of the Fermionic vacuum to the excitation of Bosonic modes, 
as well as the coupling of Bosonic and Fermionic excitations to the EM field are formulated in Ref. P8l . We have extended 
this theory to 2D chiral superconductors with anisotropic, quasi-2D Fermi surfaces, multi-band pairing and weak disorder, to 
make predictions for signatures of anisotropic chiral and helical superconductivity based on the Bosonic mode spectrum and the 
microwave response for recent theoretical models for the superconducting state of SriRuCL l39l . Here we summarize some of 
the results from the microscopic theory that reflect the coupling between Bosonic and Fermionic degrees of freedom that are 
beyond the TDGL Lagrangian dynamics for chiral ground states. 

The dynamics of the order parameter is formulated in terms of the space- and time-dependent mean-field pairing self-energy, 
which for spin-triplet fluctuations is given by 

d(p-r 7 t)= J c/ 2 pV(p,p') J ^.f K {p',e-r,t), (15) 

where F(p,p') = Vi (T/(p)T + (p , ) + T2(p)T_(p / )) is the pairing interaction in the spin-triplet Ei u channel, / K (p,e;r,/) is the 
anomalous Keldysh pair-propagator, the energy integration is over the bandwidth of attraction near the Fermi level, — £l c < e < 
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FIG. 3: Left: Anisotropy of the excitation gap. |A(p) |, based on Eqs. 14 exhibiting deep gap minima alone [110] directions in the limit £ —> 1. 
Right: Masses and Linewidths of the E ± modes resulting from gap anisotropy calculated from the Eqs. 


20 


and 


21 


at T = 0. 


+£2 ( with Cl c <C Ef, and the momentum integration is an average over the Fermi surface defined by the Fermi momentum p. 
Below we discuss the orbital dynamics for an ESP chiral ground state of the form d(p) = d A(p), with A(p) = A K, (p), with 
Y+ (p) = ( Y_Ap) ± /Yy ( p) /\/2, and the spin component of the order parameter, d, is fixed along the direction z. The orbital 
fluctuations of the Cooper pairs are represented by two complex fields. 


t/(p; r, t) = D(r, t) Y + (p) + E (r, t ) Y (p), 


(16) 


where the notation is equivalent to that in Eq. [5] of the TDGL theory. The solution to the linearized Eilenberger equations for 
the non-equilibrium pair-propagator, in terms of the time-dependent order parameter, d(p\r,t) = dd (p; r. f), and 

the coupling of charge currents to the EM field, f\ p ■ A(r ,t) leads to coupled “gap equations” for the orbital fluctuations of the 
Cooper pairs |[39l . 


t/(p;q,ai)=i Jd 2 p'V t (p,p') |-^A(p')r? , A(p') 

1 -r 


2e 

—Vp-A(q,m 
c 


7(P') + ^(P')(® 2 - 2|A(p')| 2 - V' 2 ) d(p';q, a) - X(p')A(p') 2 d'(p';q, a) 


d'(p;q,co)=* J £/ 2 p'y,(p,p') | + ^A(p') q'A^p') 
7(P / ) + ^A(p')(to 2 — 2|A(p')| 2 - r\ a ) 


2e 


V P ■ A(q,o) 


c/ , (p , ;q,a))-A(p')A*(p')^c/(p , ;q,®) > , 


(17) 


(18) 


where d'(p;r,f) = ly(v.t) A. (p) + E i: (r.t) Y (p), = v p / ■ q generates the dispersion of both Fermionic and Bosonic excita¬ 

tions, and we have expressed the gap equations in terms of Fourier modes. Note in particular that the cross-coupling terms 
between t/(p;q, ft)) and d’(p\ q. a>) are proportional to the complex amplitudes, A(p) 2 and A*(p) 2 . The Tsuneto function [(40ft . 


A(p;co,q) = ^ 


ds tanh(/3|e|/2) 
2ni y/e 2 - |A(p)| 2 


0(e 2 -|A(p)| 2 ) 


? 7 2 + 2G)(e-ft)/2) 


rj 2 -2co(e + co/2) 


(4(e -co/2) 1 - j] 2 ) (co 2 - r; 2 ) +4T7 2 |A(p)| 2 (4(e + co/2) 2 - q 2 ) (co 2 - r/ 2 ) -F 4r/ 2 1 A(p) | 2 


(19) 


determines (i) the coupling of the EM field to the Bosonic modes, (ii) the mass shifts for the Bosonic modes D± and £ , (iii) 
finite lifetimes of Bosonic modes arising from coupling to the un-bound continuum, i.e. when M > 2min[|A(p)|], and (iv) 
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coupling of the EM field to the Fermionic spectrum, including absorption of EM radiation by creation of unbound Fermionic 
quasiparticles for h<a > 2|A(p)|. Finally, the term ^y(p) is the log-divergent integral that determines the BCS instability and 
equilbrium gap function, A(p). Thus, we can regulate the divergence and eliminate the pairing interaction, Vj, in the dynamical 
equations for d( p) and d'( p) in favor of the self-consistently determined equilibrium gap function, A(p), using the identity. 


2 Id 2 pY*{p) y(p)Fv(p) = 5 Mi v for fi , v = ±. 


Energies and Lifetimes of the modes 

The Bosonic modes are obtained from the eigenvalue spectrum of the homogeneous equations, i.e. for A = 0. In the q = 0 
limit the eigen-modes are the linear combinations D ± = D(co) ±D*(—co) and E ± = E(a>) ± £'*(—©) as in the TDGL theory, 
with D representing the phase mode and D 1 the corresponding amplitude mode. For the modes with time-reversed chirality 
we obtain 


(Aoo(®)m 2 -4A 2 Aii(cu))£' + = 0, 
(Aoo(ffl) co 2 -4A 2 [Aio(®) - An(«)])£' = 0, 


( 20 ) 

( 21 ) 


where the functions A,„„ (ft>) are moments of the q = 0 Tsuneto function. For the anisotropic E| u model with basis functions given 
in Eq. 14 with p A = cos0 and p v = sin 0, we have for the anisotropic chiral ground state, A(p) = Ae ll? I{(j>). The corresponding 
moments A,„„ are then given by 


^tun (d)) 


I 2 71 


- A z 0i^A(0;®,q = O) [7(0)] 


2 n 


cos(20)] 


2m 


( 22 ) 


For the 2D chiral p-wave ground state the gap is isotropic on the Fermi circle, in which case Aiq = Aqq = 2An = A(o), with 


*(<»)= j|A| 2 


i: 


de 


tanh ^ 


\J£ 2 — |A| 2 e 2 — (®/2) : 


-© (e 2 — |A| 2 ) , 


(23) 


leading to the degenerate E ± modes with M_ i+ = M_._ = \/2A(T), in agreement with the weak-coupling limit of the TDGL 
theory; however now valid at any temperature. 

The effects of anisotropy of the pairing state on the E ± modes are computed by solving the eigenvalue equations, Eqs. 20 21 


and [22] numerically. The minimum and maximum of the anisotropic gap function are shown as a function of the anisotropy 
parameter e in Fig. [3j 2A m ; n marks the minimum in the continuum edge of un-bound Fermion pairs. The degeneracy of the E ± 
modes is resolved by the anisotropy, however the splitting of the modes for T —> 0 and relatively weak anisotropy, e < 0.4, is 
much smaller than the prediction based on the GL parameters. The splitting of the two modes is generally smaller at lower 
temperatures, and becomes strongly suppressed by the asymmetry in level repulsion between the E ± modes and the continuum 
edge at 2A mm as the higher energy mode approaches the continuum edge. 

At sufficiently large anisotropy the continuum edge of un-broken Fermion pairs at 2A m i n intercepts the excitation energy of E ± 
modes. This opens a channel for the E ± mode to dissociate into un-bound Fermion pairs, and thus leads to an intrinsic lifetime 
for the E ± mode(s), T + = h/T + , where r 1 is the width of the E ± resonance. The latter are calculated perturbatively from Eqs. 
20|21 and 22 and are shown in the right panel of Fig. [3] onsetting precisely at an anisotropy such that M- t ± = 2A m ; n . Note also 
that close to 2A m ; n the asymmetry in the level repulsion drives M _ —>• M_ + . The large asymmetry in r reflects the different 
phase space for pair dissociation of the modes governed by ImAio(cu) A> ImAn(fi)), i.e the former is an isotropic average 
over the spectrum of un-bound Fermion pairs, whereas the latter preferentially weights regions of the Fermi surface near the gap 
maximum. Thus, two key results of a self-consistent theory of coupled Boson-Fermion excitations are: (i) the mass spliting of 
the E ± mode spectrum is strongly suppressed by the asymmetry in the level repulsion from the un-bound Fermion pairs, and (ii) 
there is a large asymmetry in the lifetimes of the E ± modes that results from the different phase space available for dissociation 
into un-bound Fermion pairs by E ± Bosons. Neither of these effects could be anticipated a priori from the TDGL theory for the 
Bosonic excitations. 


Microwave Excitation of the E ! modes 

Indeed key signatures of an anisotropic chiral ground state in SioRuOa are the excitation and decay channels for the E ± 
modes. Both depend on the charge conjugation parity of the modes. Consider an EM field incident normal on a surface of 
Si' 9 Ru 04 defined by the four-fold axis of symmetry, z, and an axis lying in the x—y plane, i.e. q _L z, and with linear polarization 
also in the x — v plane, i.e. A _L z and A _L q. The EM field couples directly to the Fermionic degrees of freedom (particles and 
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FIG. 4: Power absorption spectra normalized to the high-frequency limit of the normal-state, /^(co), for T = 0, penetration depth, A/£ = 10, 
and polarization q T A along [100] directions. For this polarization direction only the E + mode is excited. Left : weak anisotropy with 
e = 0.1. A sharp E + absorption band (blue) and a broad band of dissipation from pair dissociation that is sharply peaked at 2A max (red) are 
shown. Right: strong anisotropy with e = 0.8. A sharp E + absorption band (blue) survives weak hybridization with the un-bound continuum. 
Dissipation from pair dissociation remains sharply peaked at 2A max . 


holes), generating a current (3711391 . 

l F (q,fi)) =1V/ j dr p(evp) 


1 


ft ) 2 - T ] 2 


(1 — A(p;©,q)) 


(?v A )- 


(24) 


Note that the effects of the pairing correlations on the Fermionic contribution to the charge current - the opening of a gap in the 
Fermionic spectrum and the a.c. response of the negative energy continuum (condensate) - are encoded in the Tsuneto function, 
A(p;ft),q). For T —> 0 and low frequencies, ft) < 2A mm , only the negative energy continuum (condensate) responds as an a.c. 
supercurrent - 71 /2 out of phase with the electric field - with zero dissipation. The supercurrent, and thus the self-consistently 
determined EM field, are screened by the Meissner effect and penetrate a distance of order the London penetration depth, A. 
This length scale is typically large compared to the coherence length of the superconductor, A>^ > h/Pf- hi this limit the 
EM response is dominated by the bulk excitation spectrum. For a chiral ground state the EM field also couples directly to the 
E ± Bosonic modes as shown in Eqs. 20pT The Bosonic modes also generate a charge current. 


l B (q,ffl) = ^N f J £/ 2 p(2ev p )(v p -q)A(p;co,q)) (A*(p)d(p;q,ft))+A(p)rf'(p;q,e))) . 


(25) 


Thus, the total current, 1 = Jf + Jb, can be expressed in terms of a response function, 1/ = K,y(ft),q) Ay(q, ft)), that encodes both 
Fermionic and Bosonic contributions to the a.c. surface impedance. The response function K, ; (q, ft)) determines both dissipative 
and non-dissipative contributions to the current. In particular, the microwave power absorption spectrum can be expressed as 
a sum over the modes contributing to the Joule losses of the electric field and current that penetrate into the superconductor, 
Ps(co) = — f-c T dq Re [J(g, ft)) • E* (q, ft))]. To calculate the power spectrum requires a solution of the boundary value problem 
for the incident, reflected and transmitted EM fields. This boundary value problem determines A(q, ft)) in terms of X,- ; and the 
value of the EM field in vacuum at the surface, Bq. For a detailed discussion of the boundary solution see Refs. 1581 . 

For weak anisotropy, e = 0.1, the Bosonic modes have well defined excitation energies, h(0±{ q) = M_ ,± + 4 |q| 2 /2 m -,±> 
with M_ -t « \/2A and c+ r, |vy, and carry current for finite q at frequencies below the un-bound Fermion pair continuum, 
ft,co < 2A m j n , supported by the condensate, represented by A(p;q,ft)) > 0. Resonant excitation of the E + modes occurs over 
the frequency band, M ,± < lift) < 2A m ; n , spanned by the dispersion of the E + modes. Thus, excitation of the E ± leads to an 
absorption band that is sharply peaked near threshold as shown in the left panel of Fig. |d]for e = 0.1 and E + . Note also, the 
broad band of dissipation from dissociation Cooper pairs into un-bound Fermion pairs is sharply peaked at 2A max . For strong 
anisotropy, the E ± modes broaden into resonances. However, the E 1 resonance has a narrow linewidth due to limited phase 
space for decay into Fermion pairs. Furthermore, the excitation of ground-state Cooper pairs into un-bound Fermion pairs is 
suppressed well below 2A max . Thus, even for strong gap suppression along [100] directions there remains a strong absorption 
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resonance from the weakly damped E + Bosonic mode. Observation of the E + mode, and other signatures of the E + Bosonic 
spectrum, would provide direct evidence of an anisotropic chiral ground state in Sr 2 Ru 04 . 

Summary 

Chiral superconductors which break time-reversal symmetry necessarily belong to a higher dimensional representation of 
the crystalline point group. In the cases of Sr 2 Ru 04 , 2D 3 He-A, and UPt 3 this is a two-dimensional orbital representation. 
An important consequence is that a chiral ground state supports Bosonic excitations of the time-reversed Cooper pairs. These 
excitations are degenerate for 2D 3 He-A with an excitation gap, M = \/2A, below the continuum edge of un-bound Fermion pairs. 
Crystalline anisotropy lifts the degeneracy, and for strong anisotropy can generate a low-lying Bosonic mode. Strong amplitude 
anisotropy also leads to low-lying Fermions, and thus a channel for the Bosonic modes to decay in to un-bound Fermion pairs. 
Selection rules and phase space considerations are shown to generate to large asymmetries in the lifetimes and hybridization of 
the Bosonic modes with the continuum of un-bound Fermion pairs. The excitation of the Bosonic modes by microwave radiation 
could provide clear signatures of an anisotropic chiral ground state. A detailed theory of microwave spectroscopy of anisotropic 
chiral superconductors which includes the analysis of selection rules, the effects of band structure and Fermi surface anisotropy, 
spin-orbit coupling and weak disorder will be forthcoming as a separate report ll39l . 
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